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Abstract 

The tau lepton lifetime is measured with the L3 detector at LEP using the 
complete data taken at centre-of-mass energies around the Z pole resulting in Tr = 
293.2 ± 2.0 (stat) ± 1.5 (syst) fs. The comparison of this result with the muon 
lifetime supports lepton universality of the weak charged current at the level of six 
per mille. Assuming lepton universality, the value of the strong coupling constant, 
as is found to be as{ml) = 0.319 ± 0.015 (exp) ± 0.014 (theory). 
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1 Introduction 



In the Standard Electroweak Model jlj, the couphngs of the leptonic charged and neutral 
currents to the gauge bosons are independent of the lepton generation. Measurements of the 
lifetime, r,-, and the leptonic branching fractions, B{r — > iue^r), of the tau lepton provide a test 
of this universality hypothesis for the charged current. The leptonic width of the tau lepton 

T[t £Z/£I/^) = (1) 

Tr 

2 2 5 

- (l + ep)(l + e,,,)(l + e,0, (2) 
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where I =e,fi, depends on the coupling constants of the tau lepton and the lighter lepton to 
the W boson, g^- and gi, respectively. Here rrir and mw are the masses of the tau lepton and 
the W boson. The quantities ep, ej-ad and 6^2 are small corrections resulting from phase-space, 
radiation and the W propagator, respectively. For the decay of the muon, /i Qv^i/^, the same 
formula applies, with the muon mass and coupling, and g^, replacing those of the tau. The 
comparison of the tau lifetime and leptonic branching fractions with the muon lifetime gives a 
direct measurement of the ratios gr/ ge and gr/g^i- 

Tau decays into hadrons are sensitive to the strong coupling constant, at the tau mass 
scale. Assuming universal coupling constants of the different lepton species to the W boson, 
the ratio of the hadronic width to the leptonic width, i?,-, can be expressed as: 

/2.= f^)'^-C, (3) 

where C = 1.9726 contains the same corrections mentioned above. 
Rt- is calculated in perturbative QCD [3HS]: 

ry rv rv rv 

Rr = mn<i? + \Vns?)SEw{l + - + 5.2023(^)2 + 26.366(^)3 + [l^ + d,){^Y + 5^p). (4) 

TT TT TT TT 

The quantities \Vud\ and |Vus| are elements of the Cabibbo-Kobayashi-Maskawa (CKM) quark 
mixing matrix [0], Sew [Z] and 5np IHE] describe electroweak radiative corrections and non- 
perturbative QCD contributions, respectively. The quantity is estimated as = 27.5 [Sj. 
The perturbative theoretical uncertainty is taken to be d^ = 27.5±27.5. 

This paper presents a measurement of the tau lepton lifetime with the L3 detector at LEP, 
using data taken in 1995 at the Z pole. Furthermore, data from 1994 are re-analysed using 
an improved calibration and alignment of the central tracker. The lifetime is measured from 
the decay length in three-prong tau decays and from the impact parameter of one-prong tau 
decay^ The results are combined with our previously published analyses on data collected 
from 1990 to 1993 j^. Previous measurements of the tau lepton lifetime have been reported by 
other experiments [TU] . 

2 L3 Detector 

The L3 detector is described in Ref. This measurement is based primarily on the infor- 
mation obtained from the central tracking system, which is composed of a Silicon Microvertex 

^^An A^-prong tau lepton decay indicates a decay with N charged particles in the final state. 
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Detector (SMD) ^2], a Time Expansion Chamber (TEC) and a Z-chamber. The SMD is made 
of two concentric layers of double-sided silicon detectors, placed at about 6 and 8 cm from the 
beam line. Each layer provides a two-dimensional position measurement, with a resolution of 
7 and 14 fim for normally incident tracks, in the directions perpendicular and parallel to the 
beam direction, denoted as {x, y) and z coordinates, respectively. The TEC consists of two 
coaxial cylindrical drift chambers with 12 inner and 24 outer sectors. The sensitive region is 
between 10 and 45 cm in the radial direction, with 62 anode wires having a spatial resolu- 
tion of approximately 50 /im in the plane perpendicular to the beam axis. The Z-chamber, 
which is situated just outside the TEC, provides a coordinate measurement along the beam 
axis direction. 

3 Event Sample 

For this measurement data collected in 1994 and 1995 are used, which correspond to an inte- 
grated luminosity of 49 pb~^ and 31 pb"^, respectively. 

For efficiency and background estimates, Monte Carlo events are generated using the pro- 
grams KORALZ [in] for e+e" ^ /i+/i"(7) and e_+e- r+r-(7), BHAGENE ^Mj for e+e" ^ 
e"'"e~(7), DIAG36 [12] for e"'"e~ — >• e"'"e~ff, where ff is e"*"e", r+r^ or qq, and JETSET [TH] 

for e+e~ —>■ qq(7) . The Monte Carlo events are passed through a full detector simulation based 
on the GEANT program [17^, which takes into account the effects of energy loss, multiple 
scattering, showering and small time dependent detector inefficiencies. These events are recon- 
structed with the same program used for the data. The number of Monte Carlo events in each 
process is about ten times larger than the corresponding data sample. 

Tau lepton pairs originating from Z decays are characterised by two low multiplicity, highly 
collimated jets in the detector. The selection of e+e^ t^t^{j) events is described in detail 
in Ref. [T^; here only a general outline is given. In order to have high-quality reconstruction 
of the tracks, events are accepted within a fiducial volume defined by | cos6t\ < 0.72, where 
the polar angle 9t is given by the thrust axis of the event with respect to the electron beam 
direction. The events must have at least two jets, and the number of tracks in each jet must 
be less than four. The background from e'^e" — > e+e~ (7) events is reduced by requiring the 
total energy deposited in the electromagnetic calorimeter to be less than 75% of the centre-of- 
mass energy. To reduce background from e+e~ — > (7) the sum of the absolute momenta 

measured in the muon spectrometer must be less than 70% of the centre-of-mass energy. If 
muons are not reconstructed in the muon chambers they are identified by an energy deposit in 
the calorimeters which is characteristic of a minimum ionising particle. If this is the case for 
one jet, the opposite jet is required to exhibit a hadronic signature. This rejects dimuon as well 
as cosmic-ray events. The cosmic-ray background is further reduced by requiring a scintillation 
counter hit within 5 ns of the beam crossing. In addition, the distance of closest approach to the 
interaction point measured by the muon chambers must be less than two standard deviations 
of the resolution. 

Following this procedure, 29679 and 13294 events are selected from the data collected in 
1994 and 1995, respectively. The selection efficiency in the fiducial volume is estimated to be 
76%. The purity of the tau pair sample is 98%. 
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4 Tracking 



For this measurement a high quahty of the track reconstruction is essential. A prerequisite is the 
control of the alignment between SMD and TEC and the drift time to drift distance calibration 
for the TEC. These calibrations, alignments and the estimation of resolution functions are 
performed with a clean sample of Bhabha and dimuon events, where tracks are known to 
originate from a common vertex. Particular effort is invested in the individual alignment of 
each sensor of the SMD and in the calibration of the boundaries of the TEC sectors. The 
procedure is described in detail in Ref. fH]- The performance of the track reconstruction is 
estimated from the distance between the two tracks at the vertex projected into the {x, y) 
plane. This quantity, called miss distance, is independent of the size and position of the e+e~ 
interaction region. The distributions of miss distance for Bhabha and dimuon events collected 
during 1994 and 1995 are shown in Figure A Gaussian function is fitted to both distributions, 
from which an intrinsic resolution dint = 33 /im and 31 /xm is estimated for 1994 and 1995, 
respectively. 

To guarantee good tracks for the analysis, the following cuts are made: 

• Number of hits in the TEC > 30. This ensures a good curvature measurement. 

• Number of SMD hits in the {x,y) plane > 1. This criterion selects tracks for which the 
error in the extrapolation to the vertex is well described by cxinf 

• Transverse momentum, \pt\ > 500 MeV. Tracks with lower momenta have a larger un- 
certainty in the extrapolation to the vertex due to multiple Coulomb scattering in the 



• Probability, P(x^), of the track fit larger than 1%. This requirement rejects bad fits. 



For three-prong tau decays, the decay vertex of the tau is reconstructed and its distance to the 
centre of the interaction region is measured. The decay vertex is found from a minimisation 
with respect to the vertex coordinates (x^, y^) of the following x^: 



In this equation, 5i is the distance of closest approach of a track to the decay vertex coordinates. 
The error, cx^., is the quadratic sum of the intrinsic resolution, (Tint, and the error due to multiple 
Coulomb scattering, a^asip)- Figure |21 shows the distribution of the confidence level, P{x^), of 
the fit. Decay vertices with P(x^) < 1% contain in most cases tracks with wrongly matched 
SMD hits and are rejected. The remaining distribution is flat, indicating a good description of 
the errors. 

Figure El shows the decay length distribution for the data collected in 1994 and 1995. Only 
decay length values in the range [—10, 20] mm with an error better than 5 mm are accepted 
for the lifetime determination. The number of selected decays is 4306 and 2314 for 1994 and 
1995 data, respectively. 

To obtain the average decay length, an unbinned maximum likelihood fit is applied to the 
observed decay length distribution. The likelihood function is calculated from a convolution 



SMD. 
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Decay Length Method 




(5) 
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of an exponential E, describing the tau decay time using the average decay length (L) as a 
parameter, with a Gaussian resolution function R. The likelihood function is written as: 



^ = l[i^-fB)-E®R + fB-B. (6) 

2=1 

In this equation the product runs over the accepted three-prong tau decays N^p. The second 
term on the right hand side has been added to take into account background carrying no r 
lifetime information. The background fraction, J'b, is estimated from Monte Carlo and fixed in 
the fit. The likelihood function B is evaluated from the convolution of a Dirac delta function 
with the experimental resolution function. 

The fit minimises — log£. Average decay lengths of (L) = (2.245 ± 0.037) mm and 
(L) = ( 2.265 ± 0.051 ) mm are determined for data from 1994 and 1995, respectively. The 
results of the fit are represented by the solid lines in Figure El 

The tau lifetime and average decay length are related through the following expression: 

Using this equation and taking into account the effects of radiation, lifetimes of t^- = 292.5 ± 
4.9 fs and = 295.2 ± 6.6 fs are obtained for the two data sets, where the errors are statistical 
only. 

In order to check the decay length method, the analysis is repeated on a sample of Monte 
Carlo events. The difference between the input tau lifetime and the result of the fit is assigned as 
a systematic error due to the method. Systematic effects due to the non-ideal description of the 
resolution function are estimated from a one a variation of its parameters. The deviation of the 
SMD radius from its nominal value is measured using Bhabha and dimuon events by minimising 
the track coordinate residuals from overlapping silicon sensors. Deviations of (—5 ± 5) /im 
and (—3 ± 5) yum are found in the data of 1994 and 1995, respectively. These results are 
compatible with zero, and their uncertainties are translated into a systematic error on the 
lifetime. Uncertainties in the fraction of background events carrying no lifetime information 
are estimated from a ± 50% variation of the fraction. Finally, the systematic error due to the 
fit range is estimated from the combined data by a 10% variation of their lower and upper 
bounds. This error also includes the effect of a variation of the cut on P(x^) between 0.1 and 
1.5%. Table n summarises the systematic errors for the decay length method. 

The lifetime measurements from the 6620 three-prong decays from 1994 and 1995 are com- 
bined. The systematic error due to the resolution function description is taken to be uncorre- 
lated; for the other errors a 100% correlation is assumed. The result is = 292.9 ± 3.9 (stat) ± 
1.9 (syst) fs. 



6 Impact Parameter Method 

A second measurement of the tau lifetime is obtained from the impact parameter in the plane 
perpendicular to the beam axis for one-prong tau decays. The impact parameter is the distance 
of closest approach, 6, of the track to the tau lepton production point, which is estimated by the 
beam position. In order to be sensitive to the lifetime, a sign is given to the impact parameter, 
according to the position of the intersection between the track and the tau direction of flight 
with respect to the beam position. The uncertainty on the impact parameter, (Tip, is described 
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as the quadratic sum of the intrinsic detector resolution, (Tint, the size of the interaction region, 
{cx, o'y), and a momentum dependent multiple Coulomb scattering contribution, o"ms(p), 

f^ip = o^int + ^xSin^0 + t^ycos^^ + crLip), (8) 

where is the azimuthal angle in the plane perpendicular to the beam axis. The average 
values for the interaction region size are determined from Bhabha and dimuon data and listed 
in Table El 

In contrast to the decay length method a more complicated function for the description of 
the tau decay is expected here, since for a one-prong tau decay the decay vertex is a priori not 
known. From a Monte Carlo study of the impact parameter distribution at generator level it is 
found that this function can be described in terms of three exponentials for positive and three 
exponentials for negative impact parameter values 

U{6) = {l-W)j:^e^^ +wj:^e^. (9) 

i=i i=\ \ 

In this equation, W represents the fraction of negative impact parameter values, which originate 
from an imperfect reconstruction of the tau flight direction. The slopes of the exponentials, Aj, 
contain the lifetime dependence of the distribution. 

As for the decay length method, the lifetime is extracted from an unbinned maximum like- 
lihood fit to the observed distribution. The likelihood is now determined from the convolution 
of a double Gaussian resolution function, which is obtained from Bhabha and dimuon sam- 
ples, with the function of Eqn.(9). The fit also accounts for background carrying no lifetime 
information. Figure |3] shows the impact parameter distributions from tau decays collected in 
1994 and 1995. Impact parameters with a value in the range [—0.9, 1.35] mm and an impact 
parameter error of less than 250 /xm are accepted for the measurement. The fit yields a tau 
lifetime of = 292.7 ± 3.3 fs and = 295.0 ± 4.9 fs for the two data samples. The errors 
are statistical only. 

The method is checked on a Monte Carlo sample, from which the lifetime is determined 
in the same way as for the data. The difference between the input tau lifetime and the result 
of the fit is assigned as a systematic error due to the method. Systematic effects due to the 
uncertainty of the resolution function are estimated from a variation of its parameters according 
to their errors, with correlations taken into account. The beam spot size is varied according to 
its statistical errors. The change in the central value is assigned as a systematic uncertainty. 
The effect of the average SMD radial position uncertainty is treated in the same way as in 
the decay length analysis. The systematic effect due to the knowledge of the function U{6) 
is evaluated by taking into account its statistical uncertainty and its dependence on the tau 
lifetime in the range from 250 to 350 fs. The uncertainty arising from the fraction of background 
events is estimated from a ± 50% variation of this fraction. The systematic error induced by 
the choice of the fit ranges is estimated from the combined data sample by a 10% variation of 
their bounds. Table El summarises the systematic errors for the impact parameter method. 

The lifetime measurements from the 58656 one-prong decays from 1994 and 1995 are com- 
bined. The systematic error due to the resolution function is taken to be uncorrelated. For the 
other errors a 100% correlation has been assumed. The result is r^- = 292.8 ± 2.7 (stat) ± 
2.0 (syst) fs. 
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7 Discussion 



The combination of the results obtained by the two methods with our previous ones jH] yields 



Correlations within the systematic errors are taken into account. This result supersedes all 
previous results [HIE]. This value is in good agreement with the current world average pUj . 

The measurements of the branching fractions, B{t — >■ eVeVr) = (17.806 ± 0.129)% and 
B{t —>■ fiUfj^Ur) = (17.341 ± 0.129)% [21] together with this lifetime measurement and the 
muon lifetime 120] yield Qt/Qc = 0.996 ± 0.006 and Qr/gfj. = 0.996 ± 0.006 supporting the 
universality hypothesis. 

From the tau lifetime, the tau mass, the muon mass and muon lifetime, Rr is found to be 
Rr = 3.595 ± 0.048. This corresponds to 



The first error is due to the errors of the tau lifetime and the CKM matrix elements |20j . 
The second error is the quadratic sum of the uncertainties resulting from the renormalisation 
scale, the term fourth order in as, the electroweak corrections Sew, and the non-perturbative 
correction, S^p. The renormalisation scale uncertainty is estimated following Ref. |22j by a 
variation between 0.4 < rn^/ii^ < 2.0 and is the dominant contribution to the error. Other 
contributions to the theory error as discussed in Ref. |23| are not considered. This result is in 
good agreement with other measurements of as at the tau mass jSOlElj. The value of as{m^) 
is extrapolated to the Z mass using the renormalisation group equation with the four loop 
calculation of the QCD /5-functions |2E]. The resuh, as{mz^) = 0.1185 ± 0.0019 (exp) ± 
0.0017 (theory), is in good agreement with the current world average value |20.J . 
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Table 1: Systematic errors for the decay length method 
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Table 2: The size of the interaction regions 
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Table 3: Systematic errors for the impact parameter method 
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Figure 1: Miss distance distributions from 1994 (left) and 1995 (right); Bhabha and dimuon 
events are shown in hnear scale (upper) and logarithmic scale (lower). Dots are data and the 
solid line is the result of a fit with a Gaussian. 
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Figure 2: Confidence level, -P(x^), of the secondary vertex reconstruction. 
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Figure 3: Decay length distributions from 1994 (left) and 1995 (right); Three-prong tau decays 
are shown in linear scale (upper) and logarithmic scale (lower). The hatched areas represent 
the distributions of background events carrying no lifetime information. 
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Figure 4: Impact parameter distributions from 1994 (left) and 1995 (right); One-prong tau 
decays are shown in hnear scale (upper) and logarithmic scale (lower). The hatched areas 
represent the distributions of background events which carry no lifetime information. 



16 



